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Abstract 

It is shown that the gravitational potential outside an inhomogeneous ellip- 
soid of revolution (spheroid) whose isodensity surfaces are confocal spheroids is 
identical to the gravitational potential of a homogeneous spheroid of the same 
mass. 



1 Introduction 

The result presented in this note was published in our paper entitled "Rotation 
of gas above the Galactic disk" (Astrophysics, 1988, 28, 57-64; see also Pro- 
ceedings of the Joint Varenna- Abastumani International School and Workshop 
on Plasma Astrophysics, 1986, ESA SP-251, 551-556). (Below we give Sect. 2 
of this paper, as it appeared in Astrophysics.) Recently B.P.Kondratyev (2003, 
The Potential Theory and Equilibrium Figures 1 , Moskow-Izhevsk: The Insti- 
tute of Computer Science, 624 p.) showed that the same result is valid in the 
more general case of confocal triaxial ellipsoids. 

2 The Gravitational Potential of the Disk 

In this section, we obtain the gravitational potential of an ellipsoid of revolution 
having as isodensity surfaces confocal ellipsoids. The gravitational potential of 
a body of arbitrary shape is 

where the integration is over the complete volume of the body. If the body pos- 
sesses some symmetry (exact or approximate), then the most effective method 

x We thank O.G.Chkhetiani for bringing this monography to our attention. 
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for finding the potential is an expansion with respect to orthogonal functions. 
The particular choice of the orthogonal system of functions depends on the 
symmetry. 

We shell assume that the stellar disk of the Galaxy has the form of an oblate 
ellipsoid of revolution (spheroid) . The section of the disk perpendicular to the 
plane of rotation has the form of an ellipse. It is then natural to make all the 
calculations in a system of oblate spheroidal coordinates: 
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c£v , < £ < oo, -1 < rj < 1 , < (j> < 2n . 
The Lame parameters in this system of coordinates have the form 
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h^ = c[(e + i)(i-v 2 )] l/2 - 



and the element of volume can be expressed in terms of the introduced coordi- 
nates as follows: 

d 3 r = c 3 (Z 2 +r] 2 )dZdr]d<j), (3) 

where c = V a 2 — b 2 is the half-distance between the focuses of the spheroid, and 
a and b are, respectively, the semimajor and semiminor axes of the spheroid. 
The boundary of the spheroid is determined by the relation £ = £o = b/c. 

In the oblate spheroidal coordinates we can expand R^ 1 in a series in asso- 
ciated Legendre functions [7,8]: 
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where 



e m = 1, for m = , 
e m = 2, for m > , 



P™ and Q™ are associated Legendre functions of the first and second kinds, 
respectively. 

Substituting the expansion (4) for the case £ > £o > £ , i- c -, outside the disk, 
in the general formula (1) and taking into account (3), we obtain a representation 
of the gravitational potential of the ellipsoid in the oblate spheroidal coordinates: 
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The integration over £ is from to £ , over 77 from —1 to +1, and over <fi from 
to 2tt. 

We assume that the density does not depend on <j) (axial symmetry), and 
then all integrals with m > are zero. The expression for the gravitational 
potential simplifies, 



i ]T(2n + l)P n (v)Pn(iZ)Pn(v)Qn(iO 
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We consider the terms of this series 

where $o, $1, $2, etc., are determined by 
J $ (^,r?) = -c 2 G/ arctg(l/O , 

\ Jo = 2 7 r/ ( f /! 1 p(e',r ? ')(e' 2 + '?' 2 ) J Po(r7') J Po(^K'^'; 

J = -3c 2 GJ 1?7 [£arctg(l/0 - 1] , 

\h = 2 m jfj 1 _ 1 p(^r ] ')^ 2 + V 2 )P 1 (r 1 ')P 1 ^')d^dr,'; 

\ I 2 = 2tt/ £o J 1 _ 1 p(^v)(C' 2 +v' 2 )P2(v)P2^')^'dr ] ' . 



(5) 



(6) 



(7) 



The constants Iq,I\,I2, etc., can be readily determined by specifying the par- 
ticular form of the density function. For p = i.e., when the density is 
distributed over confocal spheroids (and also in the special case p = const), all 
the constants except Jo an d h are zero. This can be readily seen by representing 
£ 2 + i] 2 in the form 



e+v 2 ^Po(v){( 
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and by using in the integration over 77 the orthogonality property of the Legendre 
functions: 

2 
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In the well-known (see, for example, [9]) special case p = const, the constants 
Jo and I 2 have the form 



h = yp£o(l + £ 2 ) = ^PV=^, 

h = -\h. 
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In (8) V = (4ir/3)a 2 b is the volume of the spheroid, and M is its mass. 
After substitution of (6), (7), and (9) in (5) we obtain 

= -^|arctg(l/0 + J(3»y 2 -l)[(l + 3e 2 )arctg(l/0-3^]} . (10) 

The oblate spheroidal coordinates £ and r\ are related to the cylindrical coordi- 
nates r and z by 




After some manipulations the expression (10) can be shown to be identical to the 
well-known expression for the gravitational potential of a homogeneous spheroid 
[9]. However, it is found that this expression is unchanged in the more general 
case p = p(£), since the relation (9) between the constants Jo and li remains 
true. Thus, the gravitational potential outside an inhomogeneous ellipsoid of 
revolution whose isodensity surfaces are confocal ellipsoids is identical to the 
gravitational potential of a homogeneous ellipsoid of revolution of the same 
mass. Similarly, the gravitational potential outside a spherically symmetric 
mass distribution does not depend on the particular distribution of the density. 

We note that the gravitational potential of an ellipsoid of revolution for 
the special case of isodensity surfaces in the form of confocal ellipsoids was 
obtained in [10], but the identity of the obtained potential and the potential of 
a homogeneous ellipsoid was not pointed out. 

< ... > 

We are grateful to A.A.Ruzmaikin and A.M.Shukurov for numerous helpful 
discussions, to M.G.Abramyan and A.M.Fridman for their interest in the work, 
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